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This paper is devoted to the study of the influence of two parallel plates on the atomic levels 
of a Hydrogen atom placed in the region between the plates. We treat two situations, namely: 
the case where both plates are infinitely permeable and the case where one of them is a perfectly 
conducting plate and the other, an infinitely permeable one. We compare our result with those 
found in literature for two parallel conducting plates. The limiting cases where the atom is near a 
conducting plate and near a permeable one are also taken. 
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■ It has been known for a long time that the consideration of boundary conditions in the radiation field imposed, 
for instance, by the presence of material plates, may alter not only the vacuum energy, as it occurs in the Casimir 
effect [1], but also the properties of atomic systems in interaction with this field. The most common examples are the 
influence of cavities on the spontaneous emission rate, or on atomic energy levels (Lamb shift modification) or even 
on the anomalous magnetic moment of the electron (g — 2 factor). In other words, we can say that the presence of 
material walls in the vicinity of atomic systems renormalizes their transition frequencies as well as the widths of their 
spectral lines. The branch of physics which is concerned with the influence of the environment of atomic systems in 
their radiative properties is usually called Cavity QED and the above examples represent only a few of them (for a 
' review see for instance ref(s) [2-5]). 

In this paper we shall investigate how the energy levels of a Hydrogen-like atom are altered when it is placed 
in a region between two parallel plates, where at least one of them is an infinitely permeable plate. These energy 
modifications are originated from the interaction between the atom and the electromagnetic vacuum fluctuations 
distorted by the presence of the plates, as pointed out by Power [6] in 1966. For the case of perfectly conducting 
plates, this problem was firstly discussed by Barton a long time ago [7] and later on by Liitken and Ravndal [8]. 
' For the particular case where only one plate is present, the interested reader can consult ref(s) [9,10]. More recently 
some generalizations were made by Barton [11,12], and Jhe and Nha [13,14]. Cavity QED between parallel dielectric 
surfaces has also been discussed in the literature [15]. 

However, although the influence of permeable plates in the spontaneous emission rate has already been considered 
in the literature [16], its influence in the atomic energy levels has not, at least as far as the authors' knowledge. 
Our purpose here is to fill this gap in the literature. For simplicity, we shall consider the following situations: (i) a 
perfectly conducting plate and an infinitely permeable one, wich we will refer to as (CP) configuration, and (ii) two 
infinitely permeable parallel plates, wich we will refer to as (PP) configuration. The former set-up, used for the first 
time by Boyer in order to compute the Casimir effect in the context of stochastic electrodynamics [17], is particularly 
interesting since it leads to a repulsive Casimir pressure [17] (see also [18,19] and Tenorio et al [20] for the thermal 
corrections to this problem). More recently, the influence of this unusual pair of plates was also considered in the 
context of the Scharnhost effect [21,22]. Regarding the latter set-up, although it leads to the same Casimir effect as 
the usual case (two conducting plates), its influence on the radiative properties of atomic systems arc different. 

In order to calculate the desired energy shifts we shall use second order perturbation theory, regularizing the relevant 
field correlations with the aid of Schwinger's imaginary time splitting, as in ref. [8]. The results are compared with the 
usual case where both plates are perfect conductors, set-up that, from now on, we shall refer to as (CC) configuration. 

From these results we shall obtain the energy shifts for an atom placed near one single conductor plate, and near 
an infinitlly permeable one. 
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II. PLATES WITH DIFFERENT NATURE (CP CONFIGURATION) 



Let us start by considering the case where the atom is placed between a perfectly conducting plate, located at 
z = 0, and an infinitely permeable one, located at z = L (as said, from now on we shall refer to this set-up as (CP) 
configuration). For this case, the corresponding boundary conditions are: 



z x E(x,y,0,t) = 
z-B(x,y,0,i) = 



z x B(x, y, L, t) = 
z • E,{x,y,L,t) = 



In the Coulomb gauge (V.A = 0) with A a = 0, we have: 

E = -A , B = VxA. 



(1) 



(2) 



It is convenient to write separately expressions for the vector potential for the transverse electric (TE) and magnetic 
(TM) modes in the following way: 



A£*(x) = V x UT(x) 
A™(x) = V x [V x u£ M (x)] 



where we defined: 



U£ B (x) = Ne z sm(kz)e 



ikT.x 



N 



U k M (x) = e z cos(kz)e lkT 

lLJ k 

and used the following notation: x — xt + zz . The wave vector is given by 
k = (kj 1 , k z ) — (fez, ky,k z ) , k z = - =: k , 



kx 1 ky IRj 

n = 0,1,2.. 



and hence, the corresponding frequencies read: 



Wk = 



k| + (n + 1/2)tt/L 



1/2 



The normalization constant = -^/2/kj.X is obtained form the condition: 

d 3 xA£*(x).A£(x) = 47r 2 <5 AV( W ( 5(kT - k^) . 



Therefore, we can write the vector potential between the plates as: 



A(x) 



E E 



d 2 k T 1 



(2tt) 2 V2^k 



a^A^(x)e-^ k * + h.c 



\=E,M n=0 1 

where the anihilation and creation operators satisfy the well known commutation relations: 

4:Tr 2 S\y5 nn <5(k T - k' T ) , 



(3) 



(4) 



(5) 



(6) 



(7) 



(8) 



(9) 



with all other commutators being zero. 

Our purpose here is to study the effect of the vacuum field fluctuations on the energy levels of an atom placed in 
a region between the plates. With this goal, we shall use perturbation theory and assume that the fields do not vary 
appreciably in the atomic scales (dipole approximation). The first non- vanishing contributions to the energy shifts 
are obtained in second order in e from: 



Ae n = e 2 

m;k,A 



(n,0|r • E(x)|m;k, A)| s 



-n c m 



(10) 
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with |n;k, A) := \n) ® |k, A), where |n) designates an atomic state with energy e„ and |k, A), a field state with one 
photon with momentum k and polarization A. It can be shown that the perturbation caused by the magnetic field 
can be neglected [8] . In the above expression, r is the electron position operator of the atom with the origin taken in 
its nucleus. 

Separating the contributions for the energy shift (10) due to the degenerate states and non-degenerate states, 
denoted respectively by |n') and \i) we write: 



Ae (i) = _ e 2 J- l(",0|r-E(x)K,k, A)| 2 (n) 



n';k,A 



Ae (2) =e 2y |(n,0|r-E(x)|l;k,A)| 2 

" ^ P_ — Fn — l,h. 



(12) 



from the degenerated atomic states Ae^ can be written as: 



Using selection rules valid for central potentials, or properties of the vacuum field, the energy contribution coming 

(13) 



M 1J = -e 2 EEiH^i-')i 2 E™^ • 

i n' k,A 

Employing Schwinger's method of imaginary time splitting we obtain (see the Appendix) 

\(0\E x \k,\)f ^ |(0|^|k,A)| 2 1^ |(0|E T |k,A)| 2 



' G + (z) 



512£ 3 tt 



E K0|^|k,A)|' = 1 

^ w k 256i 3 7r K ' y ' 



k,A 



where we defined: 



G±(z) = 0/(3, z/2L) + 0/(3, -z/2L) - ( H (3, 1/2 + z/2L)- 



-( H (3,l/2-z/2L)+(^j ±12Ch(3) (15) 

with Of and 0? being the Huruwitz and Reimann zeta functions, respectively. As a consequence, the contribution 
coming from the degenerated levels to the energy shifts are given by: 



e 2 



512^L3 ^ 



(|<nMn'}| 2 + \(n\y\n>)f)G + (z) + 2\(n\z\n')\ 2 G-(z) 



(16) 



Let us now address our attention to the contribution coming from the non-degenerate levels. For this case, we shall 
consider separately two limiting regimes, namely, where the atom is near one plate and when it is far away from the 
plates. 

Near one of the plates, it can be shown that the dominant contribution comes from >> £ n — £i [8]. Hence, 
neglecting e n — S£ and using, as before, arguments based on selection rules for atomic transitions with spherical 
symmetric potentials, or properties of the vacuum fields, it can be shown that the two contributions Ae^ and Ae^ 
(Eq's (11) and (12)) take the same form. Consequently, using the completeness of the atomic states we get: 

A £ „ = A £ W + A £ ( 2 ) 

„2 



5l2wL 3 



(n\(x A + y z )\n)G + (z) + 2{n\z z \n)G-{z) 



(17) 
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Far away from the plates (retarded regime), it can be shown that the dominant contribution comes from tu^ « 

(2) 

s n — £e [8]. Discarding now Wk, the contribution Aen becomes: 



\(n\xi 



.£|(0|E,|k,A)| 2 . 



Using the definition of the static electric polarizabilities of level n, 

IHx^)| 2 



£e - Sn 



(18) 



(19) 



as well as the matrix elements of the electric field operator obtained in the Appendix, we have in the diagonal basis 
of atomic states: 



96L 4 



. G{6) 7 \ 7 



where we defined: 



sin sin V 



(20) 



(21) 



To have the total shift away from the plates we must consider the contributions Ael 1 - 1 and Ae^ 2 ' of equations (16) 
and (20). 



III. TWO INFINITELY PERMEABLE PLATES (PP CONFIGURATION) 



Comsidering now the (PP) configuration, that is, two infinitely permeable parallel plates, the boundary conditions 
on the electromagnetic fields are now given by: 



z x B(x,y,0,t) = 
z-E(x,y,0,i) = 



z x B(x,y,L,t) = 
z • E(x,y, L,t) = 



(22) 



Using the same gauge as before (V • A = with A = 0), and writing separately the vector potential for the TE and 
TM modes, as we did for the (CP) configuration (see section II), we have: 



U k ft (x) = Ne z cos(fcz)e lkT - x 
U£ M (x) = —e z siri(kz)e lkT - x 



(23) 



mr 



k = —— , n = 0, 1, 2... k x , fc s eR 
Li 



(24) 



Following the same procedure as that employed for the (CP) configuration, we obtain after a lengthy but straight- 
forward calculation, that the energy shifts when the atom is near one of the plates are obtained from: 



64tt£ 3 



(n\(x 2 +y 2 )\n)F + (z) + 2(n\z 2 \n)F-(z) 



where we defined: 



F±(z) = Ci/(3, z/L) + Ci/(3, —z/L) ± 2C fl (3) 



(25) 



(26) 



Away from the plates, the contributions coming from the degenarated states to the energy shifts are obtained from: 

=2 



Ae« = 



647rL 3 



E 



\{n\x\n')\ 2 + \(n\y\n')\ 2 )F + (z) + 2\(n\z\n')\ 2 F.(z) 



(27) 
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and the contributions from the non degenerated states, in a basis that diagonalizes the atomic states, read: 



96L 4 



(a x +a y + a z )(F(6) + ^ j - ^a z 



where: 



F{9) 



sin 4 (#) sin 2 (0) 



(28) 



(29) 



IV. COMMENTS AND CONCLUSIONS 



This section is devoted to compare the results obtained in this paper, and the results presented in reference [8] 
where we have the (CC) configuration, that is, two conductors plates. For this boundary condition, the energy shifts 
of an atom placed near one of the plates come from: 



{n\(x 2 + y 2 )\n)F_(z) + 2(n\z 2 \n)F + (z) 



64ttL 3 

The contributions to the shifts due to the degenerated states come from: 

o2 



6471-L 3 



E 



n' _ 



n\x\n')\ 2 + \(n\y\n')\ 2 )F_(z) + 2\(n\z\n')\ 2 F + (z) 



(30) 



(31) 



for any distance from the plates. 

Far away from the plates, the enegy shifts contributions coming from the non degenerated states, in a basis that 
diagonalizes the atomic states, read: 



96L 4 



K + a y + a z ) I F(0) - ^ ) + 



(32) 



with the functions F(9) and F±(z) defined in (29) and (26) respectively. 

It is interesting to note that the function F(9) is strictlly positive along its domain (see reference [8]), but the 
function G{9) can be positive or negative, as shown in figure IV, which gives completely different behaviours for the 
energy shifts contributions (28) and (32) in comparing with (20). 

In order to do a numerical analysis of our results, let's restrict ourselves, from now on, to atoms not too highly 
excited. In this case it can be shown [8] that the contributions coming from the non degenarated atomic states, (31), 
(27) and (16), are the relevant to the energy shifts. Their signs are determined from the signs of the functions F±(z) 
and G±(z); the formers are strictly positive and the latters can change their signs, giving negative shifts for the (CC) 
configuration, while positive shifts for (PP) configuration. Note from equations (31) and (27) that these contributions 
to the energy shifts have opposite signs. Further, the roles of the longitudinal and transverse field fluctuations are 
also interchanged. For (CP) configuration, the shifts can be positive or negative. 

Here we can point out some differences between these energy shifts and the Casimir effect, another important 
manifestation of the vacuum fluctuations. For the Casimir effect, the (PP) and (CC) plates give the same attractive 
Casimir force, while for (CP) plates we have a repulsive Casimir force (but with the same L dependence as for the 
other two boundary conditions). In contrast, for the energy shifts we expect different behaviours even for those cases 
where the Casimir energies are the same, since the atom probes locally the quantum vacuum fluctuations, while the 
Casimir energy is a global quantity. 

Now we present a table with numerical results, showing the energy shifts computed for the lowest hydrogen levels 
when the atom interacts with the radiation field in vacuum state submitted to the three boundary conditions mentioned 
above. For simplicity, we assume that the atom is placed at the midle point between the plates (z = L/2). The results 
are in units of (ji(3)(am)~ 2 a/32L 3 (the results for (CC) plates can be found in reference [8]). 
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Energy Shifts 


Casimir 


double permeable 


Boyer 




-1008 


1008 





A (1) 
AAfc 210 


-576 


432 


-54 


A (1) 


-216 


288 


27 


A (1) 


-162(25+\/241) 


1296(3+\/3) 


-60,75(l+\/33) 


A (i) 


-162(25- \/241) 


1296(3-\/3) 


-60,75(l-\/33) 


A 4io 


-6156 


5184 


-364,5 


Ae (1) 


-3726 


4212 


182,2 


^ e 321 


-1782 


1620 


-60,75 


^ e 322 


-972 


1296 


121,5 



(33) 



We can see from the above table that the energy shifts for the (CC) and (PP) plates are of the same order, while 
for (CP) plates they are one order of magnitude smaller. Remark that for the Casimir effect, (CP) configuration also 
leads to a smaller force, in modulus, than (CC) original configuration, but they are of the same order of magnitude. 

As a last comment, let us analyse the limiting cases where the atom is near a unique perfectly conducting plate as 
well as a unique perfectly permeable one. For the former case we take in equation (17) the limit of the atom located 
near the conducting plate, namely, we just make z/L — > 0: 



Ae 



(n) 



647TZ 3 



(n\(x 2 + y 2 )\n) + 2(n\z 2 \n) 



(34) 



This same result can be obtained from the Lutken and Ravndal's paper [8]. 

For the second case, it would be better for calculations to have a formula that gives the energy shifts for an 
infinitely permeable plate at z = 0, and a perfectly conductor one at z = L. This expression can be obtained making 
the substitution z — > — (z — L) in equation (17): 



Ae n = - 



512tt£ 3 



(n\(x 2 + y 2 )\n)G+{-(z - L)) + 2(n\z 2 \n)G-(-(z - L)) 



(35) 



We can then take the limit z/L — > 0, giving for an atom near one infinitely permeable plate the energy shifts: 



Ae 



(«) 



647TZ 3 



(n\(x 2 +y 2 )\n) +2(n\z 2 \n) 



(36) 



We could also have obtained this expression taking the limit z/L — > in equation (25) as well. Note that the energy 
shifts for one conducting plate (34) and one perfectly permeable one (36) have opposite signs and same magnitude. 

As a last comment, we would like to emphasize the increasing importance of considering the influence of permeable 
plates in different physical situations [23,24] (see also references therein). This is clear, for instance, if we note that 
Casimir forces may become dominant at the nonometer scale and the appropriate consideration of permeable plates 
can produce repulsive forces. Recall that only attractive forces could lead to restrictive limits on the construction of 
nanodevices. 



APPENDIX 



Using the plane wave expansion of the vector potential (8), we can easily show that: 

]T|(o|£ 4 |k,A)| 2 = (o|i^|o) . (Ai) 

k,A 

These correlators are plagued with infinities and must be regularized. Choosing the Schwinger's imaginary time 
splitting, we write the regularized transverse and longitudinal vacuum fluctuations of the electric field operator 
respectivelly as: 

(0|E|(z)|0) = ]im(0|ETOz,t)E T (z,i')l ) 



6 



(0\E 2 z (z)\0) = (0|£; z (z,t)^(z,i')|0) 
1,. f / ^ 2 kT / 



cos (fc^)e 



where t' — t will be substitute by it (this is equivalent to introduce an exponential cut off). Further, defining: 

7T 



A = 



and: 



6 = \T 



l±cos(2kR) 



6 = \z 



n=0 

OO r 



l+cos( 27r(n + l/2)z/L 



n=0 



2sinh(e/2)-4 



sinh(e/2 - iff) 



c.c 



we have, omitting the limit r — > 0: 



<0|E&(*)|0) = 
(0|^W|0) = 



A 3 
4ttL 
A 3 

47TL 



<V.| • ji). )E_(r.fl) 



Expanding Eq.(A4) in powers of e, that is: 

1 e 



T±(6») + 



7 

96 V 60 



TG(9) +0(6 4 ) 



where 



T±(6) = -T-^ 



1 cos( 
6 ' sin" 



G{e)=^° s9 



sin 4 (9 sin 2 6> ' 



we obtain: 



(0|E|(z)|0) = 
<0|^ 2 M|0) = 
It's not a difficult task to show that 



vr 2 


' 8 




+ 


7 


AL 4 


?" 


T 


360 


IT 2 


' 4 






7 


4L 4 


?" 


h 2l 


360 



(0|^(z)|0) = (0|i? 2 (z)|0) = -(0|E 2 .(z)|0) 



7T 
iZ 1 



24 



G(6) 



360 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



The terms proportional to 1/e 4 in Eq's (A8) and (A9) diverge in the limit r — > 0, however they are L-independcnt, 
so that they are spurious terms with no physical significance. 

Using the same procedure adopted to compute Eq's (Al) and (A2), we can write: 
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|(0|E T |k,A)| 2 = Um y. (0|E T (z, f) |k, A)(k, A|E T (z, t')\0) 



k.A 



k.A 



1 oo 

IlimV 



n=0 



d 2 k T 1 
(2tt) 2 ^ 



Wk 



sin 2 (fcz)e^ k (*'^ 



|(0|^ik,A)i 2 _ Um y. (0\E z (z,t) \KX)(k,X\E z (z,t')\0) 



k.A 



k.A 



n=0 



/" rf 2 k T 


1 




fc 2 " 


/ (2tt)2 




w k - 


w k _ 



cos 2 (fcz)e^ k ( t '- t ) 



(A10) 



With definitions (A3) and (A4), and omitting as before the limit r — > 0, we obtain: 

^J(0|E T |k,A)| 2 _ 7T 

k.A 



CJ k 



2L 3 



1 9 

dxg-(x) - - — E_(e,0) 



E 

k,A 



l(0|^|k,A)| 2 



2L 3 



1 /*° 



(All) 



where g± (x) = {l/x i )^± (x, 9). In order to compute the above integrals we consider the analytically continued function 



9±(z,p) = -^E±{z,0) 



The integral along C p (see Fig. (2)) vanishes, which yields, using the residue theorem: 

dxg±(x,p) = (1 -e" 2 ^)" 1 - jf dzg±(z,p) + 2m ^ Resg±(x,p) 

With the definitions: 



z#0 



n=l v ' 



n=l 



and using the expansion (A6), we have: 



lim 

p^3 



(l_ e -2«pi)-i(_l) f dzg + {z,p) 
Jc. 



1 T±(0) 



3e 3 4e 



(A12) 



(A13) 



(A14) 



(A15) 



and 



\im\(l-e- 2 *n- l 2TTiY^Resg + {x,p)} 

z#0 



1 



2(4tt) 2 



-6C«(3) ± \ I C(3, 0/2tt) + C(3, -0/27r) + 



-C(3, 1/2 + 0/2tt) - C(3, 1/2 - 9/2n) + 
1 



2tt 



1287T 2 



(A16) 



which gives the desired integral: 



/oo A 00 1 

dxg±(x) = J dx-^E±(x,9) 



1 T ± (0) 1 
3^"^^ T 128^ GtW ' 



(A17) 



Subsituting this result in Eq. (All), and using the definition of S, we obtain the result (14), and a term that 
diverges in the limit t — > 0. As before it is L-independent, so that it's a spurious term with no physical significance. 
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G(0) 




FIG. 1. Graphic for G{6). 
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FIG. 2. Integration contour. 
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